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Abstract 


For  a  transverse  acceleration  wave  propagating  along  a  principal  axis  of  strain  in  a  nonlinear 
isotropic  elastic  solid,  a  simple  formula  due  to  Bicksen  relates  the  wave  speed  to  the  stress  and 
strain  state  at  the  wave  front.  We  derive  the  appropriate  generalization  of  this  result  for  finite 
deformation  viscoplasticity  models  based  on  die  multiplicative  decomposition  of  the  deformation 
gradient  into  elastic  and  plastic  parts.  The  inclusion  of  scalar  internal  state  variables  (e.g.,  to 
model  damage)  is  also  considered.  The  results  may  be  used  to  obtain  information  on  the  stress 
state  ahead  of  the  wave  if  the  strain  state  and  wave  speed  are  known.  We  discuss  applications 
to  the  analysis  of  oblique  plate  impact  tests,  where  the  transverse  wave  propagates  into  umaxially 
strained  material. 
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1  Introduction 


For  several  classes  of  isotropic  materials,  we  derive  formulas  for  the  transverse  wave 
speed  in  oblique  plate  impact  tests  in  which  a  flyer  plate  impacts  a  parallel  target 
plate  (the  specimen  under  study)  with  both  plates  inclined  relative  to  the  direction  of 
motion  of  the  flyer;  cf.  [l]-[5]-  The  velocity  imparted  to  the  target  face  upon  impact 
has  nonzero  components  in  directions  normal  and  parallel  to  the  target  face;  these  are 
referred  to  as  the  longitudinal  and  transverse  directions.  A  longitudinal  wave  brings 
the  target  into  a  state  of  compressive  uniaxial  strain  in  the  longitudinal  direction,  and 
a  slower  transverse  (or  shear)  wave  propagates  into  this  uniaxially  strained  material. 
The  problem  considered  here  is  how  to  use  the  transverse  wave  speed  to  probe  the 
stress  state  in  the  uniaxially  strained  region  ahead  of  the  transverse  wave. 

The  Cauchy  stress  tensor  is  denoted  by  T;  ai,  a2,  <73  denote  the  principal  stresses 
taken  positive  in  compression;  and  p  denotes  the  pressure: 

3 

-T  =  ^  cTiBi  O  Cj ,  p= -|trT  =  i((7i +(T2  +  cr3).  (1) 

i=l 

The  unit  eigenvectors  of  T  lie  along  the  principal  axes  of  stress.  We  assume  the 
target  plate  is  undeformed  and  stress-free  prior  to  impact;  all  measures  of  deformation 
are  taken  relative  to  this  undeformed  state.  The  principal  stretches  are  denoted  by 
Ai,  A2,  A3.  The  densities  in  the  deformed  and  undeformed  states  are  denoted  by  p 
and  po,  respectively. 

Consider  the  compressive  uniaxial  strain  generated  by  the  longitudinal  wave: 


Po/ p  —  Ai  <  1  —  A2  —  A3 .  (2) 

For  an  isotropic  material,  we  expect  that  the  principal  stresses  satisfy 

(7i  >  (72  =  (73  >  0  .  (3) 

There  is  only  one  nonzero  principal  shear  stress  (to  within  a  sign),  and  we  denote  it 
by  r: 

7"  =  |(ai  -  (72)  >  0.  (4) 

Then  by  (1),  (3),  and  (4),  the  longitudinal  stress  ctj  is  given  in  terms  of  the  pressure 
p  and  shear  stress  r  by  the  well-known  formula 

<71  =  p  +  |t  .  (5) 

Since  ai  and  Ai  can  usually  be  inferred  from  measurements  of  the  longitudinal  wave, 
we  assume  they  are  known.  Then  by  (4)  and  (5),  any  one  of  the  stresses  <72,  p,  or  r 
determines  the  other  two. 
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We  expect  that  the  transverse  wave  speed  Ut  is  influenced  by  the  stress  and 
(uniaxial)  strain  state  ahead  of  the  transverse  wave.  If  Ut  depends  on  a2,  P,  or  r 
(and  perhaps  ai  and  Ai  as  well),  then  by  measuring  Ut  and  using  the  known  values  of 
ai  and  Ai,  we  might  be  able  to  solve  for  one  (and  hence  all  three)  of  the  stresses  a2,  p, 
and  T  in  the  uniaxially  strained  region  ahead  of  the  transverse  wave.  This  approach 
is  due  to  Gupta  [2,  3];  cf.  also  Aidun  and  Gupta  [4]  and  Gonner  [5].  Our  analysis 
differs  from  that  in  [2]-[5]  in  that  we  derive  exact  formulas  for  the  transverse  wave 
speed.  The  elastic  case  is  treated  in  §2.  In  §3  we  extend  these  results  to  viscoplasticity 
models  with  instantaneous  elastic  response  based  on  a  multiplicative  decomposition 
of  the  deformation  gradient  into  elastic  and  plastic  parts.  In  §4  we  discuss  thermal 
effects  and  the  incorporation  of  scalar  internal  variables  that  affect  elastic  response; 
possible  applications  include  isotropic  damage  in  ceramics. 

Since  the  transverse  wave  may  be  structured,  the  wave  speed  Ut  of  interest  here 
is  that  of  the  wave  front,  as  this  is  the  only  part  of  the  transverse  wave  traveling  into 
a  state  of  pure  uniaxial  strain.  The  results  in  the  sequel  are  based  on  the  assumption 
that  the  transverse  wave  front  is  an  acceleration  wave.^  Across  such  a  singular  surface 
the  particle  velocity,  stress,  and  strain  are  continuous,  with  jump  discontinuities  in 
the  gradients  and  rates  of  the  transverse  component  (only)  of  these  quantities.  As 
shown  by  Abou-Sayed  and  Clifton  [1],  these  jumps  need  not  be  finite;  e.g.,  the  particle 
acceleration  immediately  behind  the  wave  front  may  be  infinite.  For  the  materials 
considered  in  this  report,  such  behavior  does  not  affect  the  formulas  for  the  wave 
speed. 


2  Isotropic  Elastic  Response 

Let  F  denote  the  deformation  gradient,  R  and  V  the  rotation  and  left  stretch 
tensors  in  the  left  polar  decomposition  of  F,  and  B  the  left  Cauchy-Green  tensor: 
F  =  VR  and  B  =  FF^  =  V^.  Then  for  isotropic  elastic  response,  T  =  T(F)  = 
'f(V)  =  r(B),  where  T  and  T  are  isotropic  functions;  cf.  [7,  §47].  The  principal 
stretches  Aj  are  the  principal  values  of  V.  The  principal  axes  of  B  and  V  are  the 
principal  axes  of  strain  in  the  deformed  state.  By  isotropy,  these  are  also  the  principal 
axes  of  stress,  so  that  by  (l)i,  V  =  Ajej  0  ej  and  B  =  ®  where 

bi  =  A?.  Let  J  denote  the  Jacobian  of  the  deformation  and  let  p  denote  the  ratio  of 

^The  other  possibility  is  that  the  transverse  wave  front  is  a  shock,  across  which  there  are  jumps 
in  the  transverse  components  of  the  particle  velocity,  stress,  and  strain.  Because  of  material  non- 
linearities,  there  may  also  be  second-order  jumps  in  the  corresponding  longitudinal  components  (cf. 
Davison  [6]  for  the  elastic  case);  for  this  reason  such  a  wave  front  is  often  called  a  quasi-transverse 
shock.  Stability  arguments  imply  that  the  shock  wave  speed  Ut  should  exceed  the  acceleration  wave 
speed  Ut-  K  the  transverse  wave  front  is  a  shock,  then  the  exact  results  in  the  sequel  may  be  con¬ 
verted  to  inequalities  by  using  the  inequality  Ut  <  Ut,  or  they  may  be  regarded  as  approximations 
if  Ut  is  replaced  with  Ut- 
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the  densities  in  the  deformed  and  undeformed  states: 


p  =  p/p(i  =  l/J,  J  =  detF  =  detV  =  \/detB  =  \Jh1h2hz  =  A1A2A3.  (6) 


Consider  a  transverse  acceleration  wave  with  direction  of  propagation  ei  and  am¬ 
plitude  (e.g.,  jump  in  acceleration)  parallel  to  62-  If  61  ^  62  (equivalently,  Ai  7^  A2), 
then  the  speed^  Ut  of  this  transverse  acceleration  wave  is  given  by  Ericksen’s  formula:^ 


<71  —  72 


71  —  72 


(7) 


where  all  quantities  are  evaluated  at  the  wave  front. 

By  (2),  (6),  and  the  isotropy  of  T,  for  compressive  uniaxial  strain  we  have 

1  <  p  =  —  1/Ai ,  A2  =  A3  =  62  =  63  =  1 5  72  =  73 .  (8) 

From  (8),  (7),  and  (4),  we  obtain  the  following  exact  formulas  for  the  speed  Ut  of  a 
transverse  acceleration  wave  propagating  in  the  longitudinal  direction  into  a  uniaxially 
strained  isotropic  elastic  material: 


pUt^  = 


(9) 


where  all  quantities  are  evaluated  at  the  wave  front.  Prom  (8)1  and  (9)  we  see  that 
7i  >  72  (equivalently  r  >  0)  is  necessary  for  the  existence  of  a  transverse  acceleration 
wave.  On  solving  (9)  for  r  and  then  using  (5)  and  (4),  we  obtain  exact  formulas  for 
the  shear  stress,  pressure,  and  lateral  stress  in  the  uniaxially  strained  region  ahead  of 
the  transverse  wave  (cf.  Scheidler  [11]): 


r=lip^-l)pUT\ 

(10) 

p  =  0-1  -  |r  =  7i  -  |(p^  -  1)pUt  , 

(11) 

7-2  =  7-1  -  2t  =  7i  -  (p^  -  1)pUt  ■ 

(12) 

For  isotropic  elastic  response,  this  solves  the  problem  posed  in  the  introduction. 

^  Ut  always  denotes  the  speed  relative  to  the  deformed  material.  For  elastic  response,  Ut  is  also 
the  ultrasonic  shear  wave  speed.  For  the  cases  covered  by  (7)  and  (9),  the  corresponding  referential 
or  Lagrangian  wave  speed  is  Ut/M- 

^Cf.  Truesdell  [8],  Truesdell  and  Noll  [7,  §74],  and  Bowen  and  Wang  [9, 10].  The  universal  relation 
(7)  does  not  require  the  existence  of  a  strain  energy  function. 
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3  Viscoplasticity  Based  on  the  Multiplicative  De¬ 
composition  of  the  Deformation  Gradient 

In  this  section  we  consider  viscoplasticity  theories  based  on  the  multiplicative 
decomposition  of  the  total  deformation  gradient  F  into  an  elastic  part  and  a 
plastic  part  F^:  F  =  F®F^;  cf.  [12,  §4]  and  [13].  F^  and  F^  have  left  and  right  polar 
decompostions  F'®  =  V'^R'®  and  F'^  =  R'^U'^,  respectively,  where  R'®  and  R'^  are 
rotations,  and  the  elastic  left  stretch  tensor  V®  and  plastic  right  stretch  tensor 
are  symmetric  positive-definite.  Then 

F  =  F^F^  =  (V^R^)  (R^U^)  =  V^RfU^ ,  R“  =  R®R^ .  (13) 

The  elastic  left  Cauchy-Green  tensor  and  the  plastic  right  Cauchy-Green  tensor 
are  defined  as  follows; 

^  (yE^2  ^  cP  ^  (F^)^F^  =  (U^)2  .  (14) 


Also  note  that  from  (13)  and  (14),  we  have 

F^  =  F(F^)-\  B^  =  F(C^)-^F^,  V^R“  =  F(U^)-^  (15) 

We  interpret  F^  as  the  deformation  gradient  from  the  undeformed  reference  config¬ 
uration  to  an  intermediate  plastically  deformed  configuration.  For  any  deformations 
for  which  F^  remains  fixed,  the  response  is  assumed  to  be  elastic,  so  that  the  stress 
T  depends  only  on  the  deformation  gradient  F®  from  the  intermediate  configuration 
to  the  current  configuration.  We  assume  that  this  elastic  response  relative  to  the 
intermediate  configuration  is  isotropic  and  unaffected  by  prior  plastic  deformation. 
Thus  there  are  isotropic  functions  T  and  T  such  that^ 

T  =  r(F^)  =  r(v^)  =  r(B®) .  (le) 

Then  T,  V^,  and  B®  are  coaxial,  so  that  by  (l)i,  =  Ef=i  Af 0  e*  and 

B®  =  Ef=i  bfei  0  Bi,  where  bf  =  (Af  )^,  and  the  principal  values  Af  of  are  the 
principal  elastic  stretches.  The  spatial  velocity  gradient  is  L  =  FF~^  (a  superposed 
dot  denotes  the  material  time  derivative),  so  by  (16)  and  (15)2, 

T  =  DT{B^)[B^] ,  =  LB^  -b  B^L^  +  F(C^)-^F^ ,  (17) 

■^By  (15)3,  is  the  left  polar  decomposition  of  Hence  if  F  and  U-^  (or  C-^) 

are  known,  then  V®  and  R#  (but  not  F^,  F^,  R^,  or  R^)  are  uniquely  determined;  cf.  Nemat- 
Nasser  [12,  §4].  Since  T  is  determined  by  V^,  it  follows  that  T  is  determined  by  F  and  (or  C^); 
this  conclusion  also  follows  from  (16)3  and  (15)2.  Thus  to  complete  the  constitutive  model,  one  need 
only  provide  evolution  equations  for  U-^  or  C^.  We  do  not  consider  specific  evolution  equations 
here,  but  instead  postulate  certain  qualitative  properties  of  these  plastic  deformation  tensors. 


4 


where  Z)T(B'®),  the  derivative  of  T  evaluated  at  is  a  fourth-order  tensor. 

We  assume  that  the  material  exhibits  instantaneous  elastic  response  in  the  sense 
that  and  are  continuous  across  an  acceleration  wave  front.  Thus 


(18) 


where  |$J  denotes  the  jump  in  $  across  the  wave  front.  These  conditions  are  equiv¬ 
alent  to  |U^]  =  [U^]  =  0  and  also  to  |(C^)-^]  =  |(C^  =  0.  Since  F  is  also 

continuous  across  an  acceleration  wave,  from  (15)2  and  (16)-(18)  it  follows  that  B® 
and  T  are  continuous  across  an  acceleration  wave,  whereas 


[t|  =  DT{B^)  [b^I  ,  |b®|  =  [L]  B^  +  B®  |Lf 


(19) 


These  relations  are  analogous  to  the  jump  relations  for  the  purely  elastic  case,  where 
=  B.  Then  an  analysis  similar  to  that  in  [7]-[10]  yields  the  following  generaliza¬ 
tion  of  Ericksen’s  formula  (7).  Assume  that  6f  7^  bf  (equivalently,  Af  7^  Af).  Then 
the  speed  Ut  of  a  transverse  acceleration  wave  with  direction  of  propagation  ei  and 
amplitude  parallel  to  62  is  given  by 


pUr^  = 


G\  —  O2 


a\  —  02 


(20) 


where  all  quantities  are  evaluated  at  the  wave  front. 

The  Jacobian  J  is  given  by  (6).  Since  F  =  F^F-^,  we  also  have  J  =  J^J^,  where 
the  elastic  Jacobian  J®  and  a  plastic  Jacobian  are  given  by 

=  det  F^  =  det  =  Af  Af  Af ,  =  det  F^  =  det  =  Af  Af  Af .  (21) 


The  principal  values  Af  of  are  the  principal  plastic  stretches.  The  pn^ipal 
distortional  elastic  stretches  Af  and  the  principal  distortional  plastic  stretches  Af  are 
defined  as  follows  (i,  j,  k  distinct): 


Af 

(jEjl/S 


(\f  AfV'* 

TP_  Af  _/AfAf\ 

UfAf)  ’ 

'  -  (JF)i/3  [xf  Af ; 

1/3 


(22) 


Then  Af  and  Af  are  independent  of  J®  and  J^,  respectively,  and  thus  are  measures 
of  elastic  and  plastic  distortion  only.  They  satisfy  the  constraints 

=  ^^:^=1.  (23) 

If  plastic  incompressibility  is  assumed,  then  =  1,  and  hence  Af  =  Af . 
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We  wish  to  apply  these  results  to  the  case  of  a  uniaxial  strain  history  ahead 
of  the  transverse  wave.  Isotropy  of  T  implies  that  (72  =  <73  if  Af  =  Af,  but  this 
latter  condition  does  not  follow  from  the  constituitive  assumptions  made  up  to  this 
point.  We  now  make  the  following  assumptions  (which  are  reasonable  for  isotropic 
materials):  for  a  uniaxial  strain  history  from  the  undeformed  state,  (equivalently, 
U^)  is  coaxial  with  F,  and  X2  =  Af.  Then® 


Ai  =  Af Af 


Af  =  A2 


E 


=  l/Af  =  l/Af 


(24) 


As  noted  above,  (24)2  implies  ^2  =  0-3.  The  relations  (22),  (24),  and  (8)  imply 


(25) 


On  substituting  (25)3,4  into  (20)  and  using  (4),  we  obtain  the  following  exact  formulas 
for  the  speed  Ut  of  a  transverse  acceleration  wave  propagating  in  the  longitudinal 
direction  into  uniaxially  strained  material:® 


pUt^  = 


(26) 


where  all  quantities  are  evaluated  at  the  wave  front.  We  may  solve  (26)  for  r ,  and 
then  (4)  and  (5)  for  p  and  The  results  are  given  by  (10)-(12)  with  replaced 

by  (^)^-  To  actually  calculate  r,  p,  and  from  these  formulas  we  need  to  know 

Af  in  addition  to  Ut,  P,  and  <7i.  On  the  other  hand,  if  cr2  is  known  (from  lateral 
stress  gauge  data  or  additional  constitutive  assumptions),  then  (26)  could  be  used  to 
calculate  Af ,  assuming  Ut,  P,  and  ai  are  known. 


®To  see  this,  note  that  for  uniaxial  strain,  F  =  V  is  symmetric  positive-definite,  and  hence 
so  is  F(U^)“\  since  F  and  are  coaxial.  V^R*  is  the  left  polar  decomposition  of  F(U^)“^ 
(cf.  (15)3),  so  the  uniqueness  of  the  polar  decompostion  implies  R*  is  the  identity  in  this  case,  since 
F(U^)“^  and  V®  are  symmetric  positive-definite.  Thus  F  =  with  all  three  tensors  coaxial 

(this  may  also  be  inferred  from  (15)2  and  (14)).  Then  Ai  =  Af  Af ,  which,  together  with  A2  =  A3  =  1 
and  Af  =  Af,  implies  (24). 

®The  Lagrangian  wave  speed  is  Ut/Xi-  For  compressive  uniaxial  strain  we  expect  that  Af  <  Ag  , 
which,  by  (25)  1,  is  equivalent  to  Af  <  1.  Then  the  denominator  in  (26)  is  positive,  and  r  >  0 
(equivalently  ai  >  (72)  is  necessary  for  the  existence  of  a  transverse  wave.  Note  that  (26)  reduces  to 
(9)  when  there  is  no  plastic  deformation,  i.e.,  when  Af  =  Af  =  1. 
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4  Discussion 


The  results  in  §2  and  §3  are  actually  valid  under  more  general  assumptions  than 
previously  described.  Suppose,  for  example,  that 

T  =  r(B;£i,...,£;^),  (27) 

where  Si, . . . ,  are  internal  state  variables  that  evolve  according  to 

4  =  /fe(T,B;ex,...,ejv)  if  5fc(T,B;£i, . . . ,£jv)  >  0 .  (28) 

Here  T,  fk,  and  gk  {k  =  are  assumed  to  be  isotropic  functions  of  their 

tensor  arguments  for  fixed  values  of  Si, . . . ,  Sjv-  Thus  the  response  is  isotropic  elastic 
for  fixed  values  of  the  internal  variables,  but  the  elastic  response  may  change  with 
changes  in  the  internal  variables,  so  that  the  overall  response  is  inelastic.  We  assume 
that  the  Sk  are  continuous  across  an  acceleration  wave.  Then  (28)i  implies  that  the 
rates  Sk  are  also  continuous  across  the  wave: 

I^fcl  =  0 ,  pfc]  =  0 .  (29) 

The  acceleration  wave  analysis  proceeds  just  as  in  the  elastic  case,  and  Ericksen’s 
formula  (7)  continues  to  hold;  cf.  Bowen  and  Wang  [10].  Hence  the  formulas  (9)-(12) 
for  transverse  waves  propagating  into  uniaxially  strained  material  remain  valid.  Of 
course,  the  values  of  the  stresses  in  these  formulas  will  now  depend  on  the  current 
values  of  the  internal  variables,  but  the  relations  between  r,  p,  <72,  <7i,  and  Ut  remain 
unchanged. 

The  viscoplasticity  model  may  be  generalized  by  replacing  (16)  with 

T  =  r(B^;£i,...,£^r).  (30) 

Here  the  evolution  equations  for  the  might  also  involve  the  plastic  stretch 
(as  well  as  any  internal  variables  characterizing  the  plastic  state  of  the  material). 
Conversely,  the  evolution  equation  for  (and  any  plastic  internal  variables)  would 
generally  involve  the  s*.  We  will  not  consider  specific  forms  for  these  equations. 
However,  if  we  assume  that  the  jump  conditions  (29)  and  (18)  are  satisfied  across  an 
acceleration  wave,  then  the  jump  conditions  (19)  continue  to  hold,  and  the  general¬ 
ization  (20)  of  Ericksen’s  formula  remains  valid.  Likewise,  (26)  continues  to  hold  if 
we  retain  the  assumption  that  for  uniaxial  strain  histories,  and  F  are  coaxial  and 

We  may  also  include  dependence  on  the  temperature  9  in  any  of  the  models  dis¬ 
cussed  here.  Assume  the  material  conducts  heat  by  Fourier’s  law,  where  the  thermal 
conductivity  may  depend  on  0,  p,  and  the  Sk-  Assuming  [0J  =  0  across  an  accelera¬ 
tion  wave,  it  can  be  shown  that  [^|  =  0  also;  cf.  [8]  and  [10].  Then  all  of  our  previous 
results  continue  to  hold. 
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We  conclude  by  emphasizing  that  in  the  derivation  of  our  results,  no  small  strain 
assumptions  have  been  made,  even  for  the  elastic  strains.  For  the  classes  of  mate¬ 
rials  considered  here,  this  actually  simplifies  the  analysis,  resulting  in  simple,  exact 
universal  relations  for  the  transverse  wave  speed. 
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